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PREFACE 

* 

The  Phoenix  Corporation  has  conducted  an  extensive  review  of  the  Depart¬ 
ment  of  Defense  Geomagnetics  Program.  While  the  entire  program  has  been  dis¬ 
cussed  in  a  separate  volume,  it  was  decided  to  address  geomagnetic  field 
models  separately  as  they  are  a  key  program  element. 

Over  the  last  decade,  the  utility  of  geomagnetic  field  models  or  the 
mathematical  representations  of  the  geomagnetic  field  has  expanded  con¬ 
siderably,  both  for  civilian  and  military  applications.  Such  an  increase 
has  been  accompanied  with  more  interest  in  models  and  demands  for  better 
accuracies  with  resultant  changes  in  techniques  of  model  computations  and 
input  data  requirements.  Thus,  it  is  also  timely  to  review  this  area  of 
geomagnetic  research.  Additionally,  such  a  detailed  overview  of  geomagnetic 
field  models  provides  a  vantage  point  for  reviewing  the  role  of  Project 
MAGNET  and  provides  assessment  of  the  general  utility  of  field  models. 
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The  study  of  the  Earth's  magnetic  field  is  one  of  the  oldest  subjects  in 
man's  quest  for  knowledge  about  his  environment.  It  has  long  been  known 
that  the  Earth's  magnetic  field  behaves  somewhat  similar  to  that  of  a  giant 
bar-magnet  lying  generally  in  a  north-south  direction,  located  near  the  Earth's 
rotational  axis.  Indeed,  such  a  description  of  Earth's  magnetic  field  was 
put  forth  by  Sir  William  Gilbert  (1540-1603)  in  his  book  De  Magnete  which 
started  the  field  of  geomagnetism.  Continuing  investigations  on  the  nature  of 
Earth's  magnetic  field,  however,  have  indicated  that  the  field  is  more  complex 
than  that  assumed  by  Gilbert.  For  example,  Gauas  (1839)  in  his  spherical  har¬ 
monic  analysis  of  the  geomagnetic  field  showed  that  an  accurate  description  of 
the  field  required  many  more  parameters  than  the  magnitude  and  orientation  of  a 
simple  bar-magnet  (dipole)  approximation.  Furthermore,  he  showed  that  onlv 
94Z  of  the  observed  field  could  be  accounted  for  by  sources  inside  the  Earth. 
Also,  there  are  irregular  spatial  variations  of  the  field  over  the  surface  of 
the  Earth,  and  irregular  temporal  variations  in  the  field  have  been  noted  since 
as  early  aa  1634.  As  a  result  of  such  studies  over  the  past  few  centuries,  it 
is  now  known  that  the  Earth’s  magnetic  field  is  composed  of  three  parts: 

a)  The  'main  field'  or  the  internal  field,  which  although  not  constant  in  time, 
varies  relatively  slowly  (usually  described  in  terms  of  years)  and  originates 
In  the  Earth's  core. 

b)  The  'external  field,'  a  relatively  small  fraction  of  the  amplitude  of  the 
main  field  which  varies  rather  rapidly,  partly  cyclically  and  partly  randomly, 
and  which  originates  outside  the  solid  Earth  in  the  ionospheric  and  magneto- 
spheric  regions. 
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c)  'Crustal  anomalous  field,'  smaller  than  the  main  field,  relatively  conatant 
in  time  and  place  and  the  resultant  of  local  magnetic  anomalies  in  the  near¬ 
surface  crust  of  the  Earth. 

As  the  magnetic  field  is  a  vector  quantity,  a  vector  sum  of  the  above 
three  component  fields  is  recorded  in  any  geomagnetic  measurement.  Of  all 
the  three  components,  however,  the  main  field  is  the  dominant  part  of  the  geo¬ 
magnetic  field.  This  field  varies  continuously  over  the  Earth's  surface,  with 
intensity  ranging  from  about  25,000  gammas  to  70,000  gammas.  Its  relative  and 
predictable  smoothness  is  what  makes  it  ideally  suited  for  applications  such 
as  navigation. 

Since  the  magnetic  field  is  a  vector  quantity,  it  is  necessary  to  measure 
three  of  the  seven  conventional  geomagnetic  elements  (Figure  1)  in  order  to  specli 
the  field  completely.  In  this  diagram,  X,  Y,  and  Z  are  the  three  mutually  ortho¬ 
gonal  components  of  the  total  field  vector  F,  with  the  X  axis  pointing  towards 
geographic  north.  H  is  projection  of  F  in  the  horizontal  plane,  and  D  and  I 
are  known  as  declination  and  inclination  respectively.  Lines  of  equal  inclina¬ 
tion,  declination,  total  field  intensity,  etc.,  when  plotted  on  maps  are  called 
' isoraagnecic  lines'  and  represent  the  variations  in  the  geomagnetic  field  over 
the  Earth's  surface.  Such  charts  for  the  F,  I  and  D  elements  of  the  field  are 
shown  in  Figures  2,  3  and  4  respectively.  Strikingly,  these  charts  show  little 
or  no  relation  to  changes  in  surface  geology  and  geography  such  as  mountain 
ranges,  submarine  ridges  and  trenches,  earthquake  belts  and  thus  it  is  indi¬ 
cated  that  the  source  of  the  field  lies  deep  within  the  Earth's  interior. 

Though,  as  has  been  noted,  the  Earth's  main  field  varies  smoothly  over  the 
surface,  it  is  far  from  permanent  in  time  as  geomagnetic  field  reversals  over 
geologic,  history  are  well  documented.  These  temporal  changes  in  the  main 
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Figure  1.  Relationship  of  geomagnetic  field  elements. 


field  are  viewed  .s  long  term  phenomena  and  are  usually  measured  over  a  period 
of  years.  Collectively  such  temporal  changes  in  the  main  field  are  called  the 
secular  variation.  From  the  recent  analysis  of  MAGSAT  data,  Langel  et  al  (1980) 
have  shown  that  the  Earth's  dipole  field  is  declining  at  a  rate  of  26  gammas/ 
year  and  if  the  present  trend  continued,  the  Earth's  field  would  decrease  to 
zero  within  1000  years.  Such  secular  variation,  however,  is  neither  constant 
over  time  nor  is  it  distributed  'miformally  over  the  earth's  surface.  In  all 
cases,  these  secular  variations  appear  to  be  regional  rather  than  worldwide. 

As  for  the  isomagnetic  maps,  the  rate  of  the  secular  variation  in  the  various 
elements  of  the  magnetic  field  can  also  be  represented  in  a  map  form  termed  iso- 
poric  charts.  An  example  of  such  a  chart,  the  secular  variation  of  the  total 
field  intensity,  F,  is  ahown  in  Figure  5. 

Though  the  geomagnetic  field  has  been  studied  for  centuries,  the  underlying 
cause  of  it  is  far  from  certain.  However,  the  most  commonly  accepted  theory  is 
that  of  a  self-exciting  dynamo,  a  concept  originally  proposed  by  Larmor  in  1919 
and  later  developed  by  Bullard  in  1949.  The  Earth's  outer  core,  which  extends 
from  a  depth  of  2900  km  to  5100  km,  is  known  to  have  the  properties  of  a  liquid 
from  seismic  evidences.  For  several  reasons,  it  is  assumed  that  this  material 
ia  a  contoination  of  iron  and  nickel,  both  good  conductors  of  electricity.  The 
self-exciting  dynamo  theory  suggests  that  the  internal  field  is  maintained  by 
electric  currents  within  the  Earth  induced  by  the  conducting  material  in  the  core 
which  is  set  into  motion  by  convection.  The  eecular  variation  of  the  main  field 
can  also  be  explained  by  this  theory  as  it  is  most  likely  connected  with  the 
changes  in  the  convection  currents  in  the  core,  in  the  core  mantle  coupling,  and 
in  the  rotational  speed  of  the  Earth. 
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The  external  part  of  the  geomagnetic  field,  the  external  field,  originates 
from  sources  outside  the  Earth  and  comprises  roughly  1Z  of  the  total  field 
observed  at  the  surface.  Ircluded  in  this  is  the  field  induced  in  the  Earth's 
crust  due  to  the  time  variability  of  the  external  field  sources.  The  variation 
in  time  of  this  field  is  much  more  rapid  than  the  secular  variation  of  the 
Earth's  main  field.  Numerous  studies  of  this  field  have  well  established  that 
these  are  associated  with  the  electric  currents  in  the  ionospheric  and  magneto- 
spheric  regions.  Though  the  external  fields  are  highly  erratic  and  range  in 
amplitude  from  several  milli-gammas  to  hundreds  of  gammas,  some  of  the  well- 
documented  effects  are  listed  below. 

i)  An  11-year  cycle  of  variations  in  the  field  intensities  are  noted 
which  correlates  well  with  the  sunspot  activities  and  has  a  latitude-dependent 
distribution. 

ii)  Diurnal  variations,  which  have  a  periodicity  of  about  a  day  and 
amplitude  averaging  about  25  gammas.  These  variations  vary  with  latitude  and 
season,  and  are  controlled  by  action  of  the  sun  on  the  ionospheric  currents. 

In  general,  two  types  of  variations  are  noted:  the  'quiet  day'  variation  (Sq) 
and  the  'disturbed  day'  variation  (S^) .  The  quiet  day  variation  is  smooth, 
regular  and  low  in  amplitude  whereas  the  disturbed  day  variation  is  less 
regular,  higher  in  amplitude  and  is  associated  with  magnetic  storms. 

iii)  Lunar  diurnal  variation  (L)  with  a  periodicity  of  approximately 
25  hours  and  having  amplitudes  about  one-f if taenth  the  amplitude  of  the  solar 
Sq  diurnal  variation.  These  vary  cyclically  through  the  month  and  are  asso¬ 
ciated  with  lunar  gravitational  influence  on  the  ionosphere. 

iv)  Magnetic  storms  which  are  transient  disturbances  lasting  for  several 
days  and  having  amplitudes  as  large  as  1000  gammas  in  most  latitudes  and  even 
greater  in  polar  regions,  where  they  are  usually  accompanied  by  aurora.  Magnetic 


storms  are  usually  not  predictable,  but  they  tend  to  come  at  interval*  ot 
27  days,  a  frequency  correlating  with  the  rotation  of  the  sun. 

A  more  detailed  breakdown  of  the  general  characteristics  of  external 
magnetic  fields  is  given  in  Table  1  which  lists  the  various  types  of  external 
field  phenomena  observed  over  the  Earth's  surface,  its  typical  period  and 
amplitude  and  its  probable  origin. 

The  third  component  of  the  geomagnetic  field  is  the  crustal  anomalous 
field  which  arises  as  a  result  of  the  variations  in  the  magnetic  mineral 
content  of  crustal  rocks.  Such  fields  are  essentially  time  invariant  and 
local  in  nature  and  are  also  called  local  magnetic  anomalies.  The  magnitude 
of  these  anomalies  is  typically  on  the  order  of  several  hundred  gammas  and 
may  extend  from  one  to  several  hundred  square  kilometers  in  areal  coverage. 

Such  magnetic  anomalies  are  the  primary  targets  of  exploration  geophysical 
magnetic  surveys.  The  depths  to  the  sources  of  these  anomalies  are  presumed 
to  be  less  than  20  km,  as  che  Curie  Point  temperature  is  exceeded  below  this 
depth.  Thus,  the  crustal  anomalous  fields  are  primarily  associated  with  the 
near-surface  features. 

It  can  thus  be  summarized  that  the  Earth's  magnetic  field  characteristically 
arises  from  two  classes  of  sources,  one  internal  to  the  Earth  and  the  other 
external  to  it.  The  part  internal  to  the  Earth,  is  a  relatively  slowly  time 
varying  field  with  time  constants  of  tens  to  thousands  of  years  and  can  thus 
be  regarded  as  a  permanent  or  steady  field  for  many  applications.  This  steady 
part  of  the  Earth's  total  field  is  also  known  as  the  main  field  and  it  results 
primarily  from  convective  motion  in  the  core.  Approximately  90X  of  the  main 
field  is  dipolar  in  nature,  with  the  ideal  dipole  centered  close  to  the  geo- 
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graphic  center  of  the  Earth  and  its  axis,  inclined  about  11.5  degrees  to  the 
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GENERAL  CHARACTERISTICS  OF  EXTERNAL  MAGNETIC  FIELD  TEMPORAL  VARIATIONS 
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axis  of  rotation.  The  remaining  10Z  of  the  main  field,  often  termed  the  residual 
field,  is  nondipolar.  It  consists  of  both  large-scale  anomalies  (up  to  thousands 
of  kilometers)  believed  to  be  generated  by  eddy  currents  in  the  fluid  core,  and 
small-scale  irregularities  (up  to  tens  of  kilometers)  originating  from  residual 
or  induced  magnetism  in  the  Earth's  crust.  The  part  of  the  Earth' 8  magnetic 
field,  originating  from  outside  the  Earth,  1b  called  the  external  field  and 
recent  analyses  have  indicated  it  to  be  less  than  1Z  of  the  total  field  measured 
at  the  Earth's  surface. 

Historically,  the  use  of  the  Earth’s  magnetic  field  has  been  in  the  field 
of  navigation  relying  on  the  directive  action  of  the  field  upon  compass  needles. 
For  such  purposes,  declination  maps  were  prepared  by  hand  by  drawing  smooth 
curves  through  the  measured  data  available  from  permanent  observatory  locations 
and  periodic  ship  measurements.  However,  as  the  knowledge  about  the  Earth's 
field  and  its  applications  in  other  areas  of  physics  widened,  it  became  Increas¬ 
ingly  apparent  that  such  a  chart  or  map  representation  of  the  field  was  not  ade¬ 
quate.  For  example,  in  space  sciences  applications,  it  is  necessary  to  know 
the  Earth's  field  at  satellite  altitudes  and  the  charts  which  represent  field 
values  measured  at  the  Earth's  surface  are  of  little  use.  Also,  in  magnetic 
surveys  for  mineral  or  petroleum  exploration,  it  is  necessary  to  remove  the 
background  or  the  regional  field,  so  that  the  anomalous  field  associated  with 
the  geological  structures  can  be  detailed.  Thus,  a  mathematical  representation 
of  the  Earth's  field  suitable  for  such  applications  is  required.  Such  represen¬ 
tations  not  only  help  extrapolate  field  values  at  locations  where  specific 
measurements  have  not  been  made,  but  they  are  also  useful  in  assessing  the 
possible  mechanisms  or  sources  of  the  main  field. 
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Such  a  description  and  representation  of  the  geomagnetic  main  field  by  a 
mathematical  model  is  called  a  geomagnetic  field  model.  Such  field  models  are 
usually  derived  by  fitting  a  spherical  harmonic  series  to  a  set  of  global  data, 
a  method  first  put  forth  by  Gauss  in  1839.  In  such  a  spherical  harmonic  series 
representation  of  the  field,  usually  the  terms  pertaining  only  to  the  sources 
interior  to  the  Earth  are  considered,  as  the  contributions  from  external  sources 
are  less  than  li  and  are  lost  in  the  overall  accuracy  of  the  fit  of  the  model. 
Thus,  these  models  are  also  referred  to  as  internal  field  models. 

With  the  advent  of  satellite  measurements,  quantitative  modeling  of 
external  fields  has  also  been  attempted.  Hovever,  such  efforts  have  not  been 
very  successful  primarily  because  of  the  difficulty  in  separating  the  internal 
and  external  fields  from  the  measured  data  due  to  their  relatively  small  ampli¬ 
tudes  and  extreme  temporal  and  spatial  variations.  Nonetheless,  external  field 
models  are  becoming  of  increasing  importance,  as  the  accuracies  of  magnetic 
surveys  have  reached  a  point  where  a  more  accurate  regional-residual  separation 
in  crustaJ  anomalies  studies  is  desired.  Furthermore,  since  satellites  are  row 
being  used  to  map  the  larger  crustal  anomalies,  models  of  external  fields  are 
needed,  as  at  satellite  altitudes  the  external  field  effects  are  much  more 
signif leant. 

The  use  of  such  field  models  in  many  different  areas  of  geomagnetic  studies 
is  manyfold  and  many  such  models,  with  different  data  bases,  of  different  time 
validity  and  of  different  order  and  degree  of  spherical  harmonic  series 
expansion  have  been  published  in  the  literature  since  Gauss's  first  model  in 
1839.  Many  excellent  discussions  on  field  models  themselves  and  comparative 
studies  of  various  field  models  have  also  appeared  in  the  literature  in  the 
recent  years  (Coles,  1979;  Mead,  1979;  Barraclough,  1978;  Barraclough  et  al., 
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1978;  Earraclough,  1976;  Ragan  and  Cain,  1975;  Kane,  1973;  and  Cain,  1971). 
This  report  builds  on  sur.h  studies  and  presents  a  brief  summary  of  the  methods 
of  determining  field  models  end  the  current  atate-of-the-art  in  field  model 
studies. 
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FIELD  MODELS 


Definition 

A  field  model  is  basically  a  four  dimensional  function  representing  the 
geomagnetic  field  at  any  point  in  space  and  time.  It  is  a  function  of  colati¬ 
tude  4>,  longitude  G,  geocentric  distance  r  and  time  t.  Since  the  geomagnetic 
field  is  derivable  from  a  potential  function,  and  because  any  potential  could 
be  expressed  in  terms  of  spherical  harmonic  series  expansion,  the  field  models 
themselves  are  usually  expressed  in  terms  of  spherical  harmonic  aeries. 

The  scalar  geomagnetic  potential,  V,  at  an  external  point  (r,  0,  $)  arising 
from  sources  inside  the  earth  could  be  expressed  in  terms  of  spherical  harmonic 
series  a9  (e.g.,  see  Chapman  and  Bartels,  1940): 

n 

53  (g0  cos  m4  +  sin  mb)PB(0),  (1) 

m=0  n 

where 

a  «  mean  radius  of  the  earth; 

m,n  *  order,  degree; 

pm(0)  -  Schmidt's  quasinoraalized  spherical  functions;  and 

n 

gm,  hm  -  Gauss  coefficients 

n  n 

The  three  orthogonal  components  of  the  field  could  be  derived  by  taking 
the  negative  gradient  of  the  potential.  The  northward,  eastward,  and  downward 
components  of  the  field  are  thus 

X  -  (1/r)  (3V/30) 

Y  -  -(1/r  sin  0)  (3V/3$)  (2) 

Z  -  3V/3r 
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Carrying  out  the  differentiation  of  equation  (1),  the  X,  Y ,  Z  components 
are  given  explicitly  as 

n+2 


x  -  E 


*  n  /  \  i 


(gm  coa  n4  +  hn  sin  m$)  ~  pm(Q) 
n  n  n 


■  k  s© 


n+2 


(gm  sin  -  hm  coa  m$)  p°0) 

n  n  n 


(3) 


k  k  @ 


n+2 


(n+1)  (g™  coa  +  h™  sin  m$)  p®(0) 


It  is  thuB  clear  from  the  above  expressions  that  if  the  set  of  Gauss 
coefficients  are  known  in  the  series,  the  field  components  or  any  other  element 
of  the  field  may  easily  be  computed  at  any  point  over  the  earth's  surface.  Such  a 
set  of  Gauss  coefficients  are  normally  derived  or  adjusted  by  a  least  squares 
procedure,  wherein  some  quantity  such  as  the  weighted  sum  of  the  squares  of  the 
difference  between  the  value  calculated  from  the  model  and  the  observed  data  is 
minimized;  the  sum  being  taken  over  a  network  of  measurements  covering  the  entire 
globe.  It  is  to  be  noted  here  that  if  the  global  data  set  consists  entirely  of 
field  components  only,  the  set  of  Gauss  coefficients  can  be  derived  hy  the  least 
squares  method  in  a  straight-forward  manner.  If,  however,  the  date  cot  also 
consists  of  some  other  elements  of  the  field,  for  example,  scalar  total  field, 

F,  then  a  non-linear  least  squares  procedures  must  be  used,  as  the  resulting 
normal  equations  are  non-linear  in  the  Gauss  coefficients.  In  such  a  case,  the 
equations  are  quasi-normalized,  and  an  assumed  initial  estimate  of  the  aet  of 
coefficients  are  adjusted  iteratively,  until  a  fit  to  desired  accuracy  is 
achieved.  Normally,  this  latter  procedure  is  used  in  the  field  model  calcula¬ 
tions,  since  the  global  data  set  Is  invariably  composed  of  other  elements  as 
well  as  the  field  components  themaelvt  s. 
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In  the  spherical  harmonic  series  representation  of  the  geomagnetic 
potential  (expression  (1)),  it  is  noted  that  the  order,  n  of  expansion  should 
be  carried  out  to  infinity  for  a  perfect  theoretical  representation.  In 
practice,  however,  the  expansion  to  such  a  high  degree  is  neither  possible 
nor  warranted.  Each  value  of  n  in  the  spherical  harmonic  expansion  is  a 
global  wave  number,  i.e.,  etch  harmonic  represents  variations  of  the  poten¬ 
tial  whose  wavelengths  are  approximately  40,000/n  km.  Figure  6  shows  this 
relationship  graphically.  It  can  be  noted  from  this  figure,  that  by  carrying 
out  the  expansion  to  degrees  higher  than  8  to  10,  the  representation  of  smaller 
wavelength  features  in  the  model  improves  only  very  slowly.  Also,  as  the 

degree  n  increases  in  the  model,  the  number  of  coefficients  increases  rapidly 

2 

as  [(n-t-1)  -  1],  Hence,  the  number  of  calculations  and  the  amount  of 

associated  computer  time  ^nd  storage  used  in  the  least-squares  analysis 
increases  greatly,  approximately  as  the  fourth  power  of  n.  Thus  the  series 
is  truncated  usually  somewhere  between  the  harmonic  degree  of  8  to  12,  depend¬ 
ing  upon  the  accuracy  of  the  data.  Studies  have  indicated  that  field  models 
of  such  degree  and  order  are  sufficient  to  model  the  field  arising  from 
the  sources  in  the  core  of  the  earth.  Some  field  models  of  even  higher 
degree  and  order,  however,  have  also  been  proposed. 

A  time  parameter  for  modeling  secular  change  can  be  readily  introduced 
into  the  model  by  expanding  the  Gauss  coefficients  in  a  finite  Taylor  series 
about  some  mean  time  of  the  data  set  termed  the  epoch  tQ: 
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Figure  6.  Appropriate  relationship  between  the  d 
harmonic  series  expansion  and  the  wave 
represented  in  the  model. 


g 


"  8C  +  g'(t  ~  tc>  (t  -  tQ)2  +  ... 

and  (A) 

h  -  hc  +  h'(t  -  tQ)  +  ^  (t  -  to)2  +  ... 

Thus,  the  least  squares  fit  is  usually  made  in  four  dimensions  (r»  6.  $>  t)  and 
the  complete  set  of  coefficients  consists  of  the  Gauss  coefficients  (g,  h)  and 
at  least  the  first  order  secular  change  coefficients  (g',  h').  However,  some 
models  are  derived  by  first  independently  estimating  the  secular  change  terms 
from  the  repeat  stations  or  observatory  data  and  then  reducing  the  entire  data 
set  to  a  common  epoch  before  least  squares  fit  is  made. 

Practical  Problems 

It  must  be  recognized  that  field  models  are  essentially  extrapolation 
functions,  i.e.,  based  on  a  set  of  observed  data,  they  are  used  in  predicting 
field  values  at  locations  over  the  global  area  where  specific  measurements 
have  not  been  made.  For  such  a  prediction  or  extrapolation  to  be  accurate 
the  nature  and  limitations  of  the  data  base  from  which  the  model  is  derived 
mu6t  be  fully  examined  and  considered. 

The  first  problem  in  accurate  model  determination  comes  from  the  fact  that 
the  available  observed  data  are  limited  and  irregularly  distributed  over  the 
global  area.  Ideally,  what  would  be  desired  is  a  uniform  and  dense  coverage 
of  data  all  over  the  globe.  However,  examination  of  charts  of  presently  avail¬ 
able  land,  sea  and  airborne  magnetic  survey  coverage  indicate  that  dense  cover¬ 
age-  is  available  in  most  parts  of  the  globe,  but  there  are  still  sizeable 
areas  where  the  coverage  Is  sparse  or  completely  nil.  Such  irregular  distri¬ 
bution  of  data  prohibits  the  accurate  determination  of  higher  order  harmonic 
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coefficients  In  the  modal  and  nay  cairte  an  aliasing  effect  In  the  lower  order 
coefficients.  It  is  one  of  the  fundamental  problems  in  spherical  harmonic  series 
field  models  that  for  calculation  of  each  harmonic  coefficient,  knowledge  of 
the  field  over  the  whole  surface  is  required  and  little  extra  accuracy  is  gained 
locally  by  having  a  very  dense  local  net  of  observations. 

In  relation  to  the  effect  of  spatial  distribution  on  the  model  calculation, 
the  Becular  variation  of  the  data  or  the  temporal  distribution  of  the  data  alao 
must  be  considered.  To  obtain  a  sufficient  amount  of  data  distributed  reason¬ 
ably  well  over  the  global  area,  measurements  obtained  over  a  considerable  peri¬ 
od  of  time  must  be  utilized.  During  such  a  period,  the  field  values  would  have 
changed  in  a  nonuniform  and  unpredictable  manner  owing  to  the  secular  variation 
of  the  main  field.  Thus,  the  model  calculation  would  be  in  serious  error.  One 
way  to  compensate  such  error  would  be  to  reduce  all  the  data  to  a  common  epoch 
by  estimating  the  secular  changes  which  occurred  during  the  time  interval  of 
data.  Such  correct  estimation,  however,  is  made  difficult  as  the  magnetic  ob¬ 
servatories,  where  the  secular  changes  are  continuously  recorded,  are  sparsely 
distributed  over  the  globe  and  moreover,  they  are  usually  restricted  to  land 
masses  leaving  the  vast  ocean  areas  of  the  globe  with  no  record  of  secular  varia¬ 
tions.  The  other  way  to  compensate  the  effects  of  secular  variation  in  the  data 
is  to  incorporate  the  time  term  in  the  field  model  S6  previously  mentioned,  by 
expanding  the  Gauss  coefficients  in  a  finite  Taylor  series  about  some  mean  time 
in  the  data  set  and  simultaneously  solving  for  the  time  derivatives,  called  the 
secular  change  coefficients.  In  recent  field  models,  this  latter  approach  it 
preferred  as  the  models  then  can  not  only  be  used  for  spatial  predictions  but 
for  temporal  predictions  as  well. 


The  problem  of  nommiform  global  coverage  was  the  driving  force  behind  the 
Initiation  of  Project  MAGNET  program  by  the  U.S.  Naval  Oceanographic  Office. 

This  is  an  ongoing  project  to  collect  high  altitude  vector  measurements  over 
most  of  the  free  world.  Although  limited  by  inherent  logistics  of  airborne 
operation  this  program  has  provided  an  invaluable  data  set  for  field  model  cal¬ 
culation.  Two  limitations  of  such  a  program,  however,  are  the  fact  that  it 
cannot  obtain  global  coverage  because  of  the  political  realities  and  that  such 
coverage  cannot  be  obtained  in  a  short  enough  interval  that  the  secular  change 
remains  constant.  This  problem  has  been  somewhat  minimized  in  the  past  decade 
with  the  aid  of  satellite  measurements.  Until  recently,  such  data  have  only 
been  of  the  scalar  total  field  and  aa  Stem  and  Bredekamp  (1974)  have  noted, 
any  model  based  on  such  data  is  very  precise  in  representing  the  total  field 
values  at  satellite  altitudes  but  may  be  deficient  in  determining  accurate 
vector  values.  The  recent  launch  of  MAGSAT  satellite  provides  scalar  as  well 
as  vector  measurements  over  the  globe  at  altitudes  lower  than  earlier  satel¬ 
lites  and  thus  these  data  are  expected  to  yield  a  better  model  for  the  current 
epoch.  However,  it  must  be  noted  that  because  of  the  lower  resolution  of 
satellite  data,  it  does  not  replace  the  measurements  from  programs  such  as 
Project  MAGNET.  The  two  data  sets  are  quite  complementary. 

Although  the  method  of  determining  the  coefficients  on  the  field  models 
is  mathematically  quite  straightforward,  the  process  of  selecting,  reducing 
and  weighting  of  the  date  is  not.  Usually,  the  data  base  for  the  model  calcula¬ 
tion  comprises  diverse  data  sets,  i.e.,  it  may  include  data  from  airborne, 
land,  marine,  and  satellite  surveys  ss  well  as  observatories.  As  each  of  the 
different  data  sets  are  subject  to  differsnt  kind  of  errors  and  have  different 
resolution,  their  accuracies  must  be  carefully  estimated  and  the  data  carefully 
screened  for  spuriouB  values.  Usually,  the  data  from  observatories  ars  given 


the  highest  weight  because  of  their  extreme  accuracies.  Such  weighting  and 
selection  of  the  data  ultimately  affects  the  utility  of  a  field  model. 

In  addition  to  such  practical  concerns,  the  basic  limitations  of  spherical 
harmonic  series  in  representing  the  observed  geomagnetic  field  accurately,  must 
also  be  considered.  As  has  been  noted,  the  Earth's  main  field  is  composed  of 
long  wavelength  anomalies  arising  from  the  core  sources  and  small  wavelength 
anomalies  arising  from  sources  in  the  crust.  These  long  wavelength  anomalies 
are  modeled  adequately  by  a  spherical  harmonic  series  expansion  of  degree  and 
order  ten,  however,  the  small  wavelength  anomalies,  which  may  be  on  the  order 
of  10-i00  km,  are  not.  To  model  anomalies,  whose  wavelength  are  on  the  order 
of  10-100  km,  would  require  expansions  up  to  degree  of  400  to  4000.  The  time 
and  storage  required  for  carrying  out  the  spherical  harmonic  series  to  such  a 
high  degree  would  be  prohibitively  large,  ruling  out  the  use  of  spherical  har¬ 
monic  series  to  include  crustal  anomalies  in  its  main  field  representation. 

Though  this  limitation  of  spherical  harmonic  series  is  not  important  in  geophysi¬ 
cal  exploration  applications,  its  impact  for  global  studies  and  apace  science 
application  is  quite  pronounced,  since  here  a  best  representation  of  the  total 
field  is  desired. 

State-of-the-Art 

The  history  of  calculation  and  representation  of  Earth's  main  field  by 
spherical  harmonic  series  goes  back  more  than  a  century  since  Guass  first 
showed  in  1839  that  such  a  field  model  could  conveniently  be  used  to  represent 
the  field  distribution  over  the  Earth's  surface.  Since  then  many  such  field 
models  have  been  proposed  and  published  in  the  literature.  Barraclough  (1978) 
has  recently  published  an  exhaustive  survey  of  over  264  field  models  that  have 


appeared  in  the  literature  up  to  1973  giving  a  detailed  account  of  the  method* 
and  the  data  bases  used  in  the  different  model  derivation* along  with  a  complete 
list  of  the  Gauss  coefficients. 

It  must  be  recognized  chat  the  calculation  of  field  models  is  an  ongoing 
process,  i.e.,  field  models  must  be  constantly  updated,  bb  more  new  data  accumu¬ 
lates  and  older  data  become  obsolete  because  of  the  unpredictive  nature  of  the 
secular  variation  of  the  main  field.  Though  the  methodology  of  representing 
the  main  field  by  spherical  harmonic  series  has  remained  the  same  over  the 
years,  two  important  advances  in  the  method  of  calculation  of  field  models 
should  be  noted.  In  the  past,  befoie  the  availability  of  modem  digital  com¬ 
puters,  the  coefficients  of  the  harmonic  series  had  to  be  computed  by  hand  using 
numerical  quadratures.  This  required  availability  of  data  distributed  in  a 
regular  fashior  in  terms  of  latitude  and  longitude.  Thus,  the  available  survey 
and  observatory  data  had  to  be  plotted  in  a  chart  or  map  form  and  had  to  be  inter¬ 
polated  to  obtain  a  regular  grid  of  data.  Considerable  smoothing  and  modifica¬ 
tion  of  original  data  was  thus  involved  resulting  in  less  accurate  field  models. 
Modern  analysis,  however,  takes  advantage  of  the  power  of  digital  computers  and 
works  directly  on  the  set  of  irregularly  distributed  data,  thus  eliminating  the 
intermediate  stage  of  grid  or  map  preparation.  The  other  advancement  in  field 
model  derivation,  again  with  the  help  of  digital  computers,  is  the  consideration 
of  the  oblateness  of  the  earth  which  is  important  for  the  true  representation  of 
the  field.  Also,  the  computational  capabilities  of  such  modem  computers  facili¬ 
tate  the  use  of  a  larger  data  base  and  expansion  of  spherical  harmonic  series 
to  a  larger  degree  and  order.  Such  improvements  in  field  model  definition  have 
resulted  in  a  better  and  more  accurate  representation  of  the  earth's  main  field 
globally.  For  example,  the  best  recent  models  are  capable  of  predicting  the 


average  quiet-time  field  over  most  of  the  Earth's  surface  with  an  accuracy  of 
about  .2  degrees  in  direction  and  200  gammas  in  magnitude  whereas  the  older 
models  could  give  at  best  one  to  two  degrees  in  the  directional  accuracy. 

Another  advancement  in  the  method  of  representation  of  the  main  field  by 
spherical  harmonic  series  has  been  in  the  terms  of  analyzing  the  secular  changes 
of  the  main  field  with  secular  acceleration  terms  now  also  being  considered 
(Barraclough  and  Malin,  1979;  Fougere,  1969;  and  Cain  et  ai.,  1967). 

As  previously  mentioned,  field  models  are  usually  computed  to  order  and 
degree  of  about  ten.  This  has  been  extended  somewhat  in  recent  years  to  twelve 
or  thirteen  as  spectral  studies  (e.g.,  Cain,  1975)  indicate  that  modeling  of 
the  core  field  requires  such  resolutio  .  Even  higher  order  models  are  desirable 
for  some  applications  but  the  computer  storage  and  time  limitations  associated 
with  spherical  harmonic  analysis  inhibit  this.  A  solution  to  the  excessive  com¬ 
puter  time  and  storage  requirement  of  the  conventional  harmonic  series  analysis 
has  recently  been  offered  by  T1.  ■>  Analytical  Sciences  Corporation  (TASC)  group. 
Like  the  Fast  Fourier  Transform  (FFT) ,  they  have  proposed  a  fast  computational 
scheme  for  the  spherical  harmonic  expansion  which  seems  to  drastically  reduce 
the  time  and  storage  requirements.  The  method,  however,  has  a  drawback  in  that 
a  gridded  set  of  data  is  required  for  the  model  construction,  and  only  gridded 
data  are  computed  from  the  model.  In  essence,  the  method  involves  expansion  of 
the  surface  spherical  harmonic  in  a  Fourier  series,  which  may  then  be  evaluated 
using  standard  FFT  technique. 

Another  approach  to  avoid  the  limited  capability  of  present  field  models 
has  been  offered  by  Brown  (1976) .  Unlike  using  the  spherical  harmonic  repre¬ 
sentation  as  is  conventionally  done  in  field  models.  Brown  has  suggested  use  of 
Walsh  functions  to  represent  the  data.  These  functions  are  binary  and  ortho¬ 
gonal,  with  a  frequency  characteristic  similar  to  that  of  standard  spherical 
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RECENT  FIELD  MODELS 


As  has  been  noted,  geomagnetic  field  modeling  is  an  ongoing  and  active 
area.  Table  2  presents  characteristics  of  some  of  che  representative  field 
models  published  in  recent  years  and  forms  the  basis  for  the  following  discus¬ 
sion.  For  details  on  earlier  models,  readers  are  referred  to  Barraclough ' s 
(1978)  paper,  where  an  exhaustive  summary  and  listing  of  models  published 
through  1973  are  presented. 

Of  all  the  models  listed  in  Table  2,  the  International  Geomagnetic  Reference 
Field  (IGRF)  model  is  probably  the  best  known  and  most  widely  used  model.  It 
emerged  primarily  as  a  result  of  an  international  need  to  have  a  standard  field 
model  to  which  all  the  data  should  be  referred  or  reduced  to,  thereby  facili¬ 
tating  inter-comparison  among  different  data  setB.  The  first  such  model  (1GRP65) 
was  chosen  at  a  symposium  held  by  the  International  Association  of  Geomagnetism 
and  Aeronomy  (IAGA)  in  Washington,  D.C.  in  1968.  The  model  was  to  represent 
the  field  at  epoch  1965.0  with  secular  change  coefficients  that  would  extend  its 
usefulness  for  at  least  a  decade.  Various  groups  of  users  were  asked  to  present 
field  models  for  consideration  by  the  IAGA  committee  and  the  IGRF65  was  developed 
as  a  compromise  best  model.  Its  main  field  coefficients  are  a  weighted  combina¬ 
tion  of  four  field  models  and  the  time  terms  ere  an  average  of  five  models.  A 
detailed  description  of  IGRF65  has  been  given  by  Cain  and  Cain  (1971). 

During  the  years  of  its  use,  however,  the  inherent  weakness  of  the  secular 
coefficient  terms  of  the  IGRF65  became  rather  evident,  i.e.,  the  predicted  value 
from  the  model  at  a  later  time  did  not  accurately  represent  the  then  measured 
field  at  many  regions  of  the  globe  (e.g.,  Regan  and  Cain,  1975;  Petkovic  and 
Whitworth,  1975).  Accordingly,  resolutions  passed  at  the  IAGA  General  Assemblies 
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in  Kyoto  1973  and  Grenoble  1975  resulted  in  the  production  of  a  revised  model 
IGRF75,  effective  1975.0.  The  IGRF75  retains  the  main  field  terms  of  IGRF65 
updated  to  1975.0  by  means  of  the  IGRF65  secular  terms.  Its  new  secular  terms 
are  the  averages  of  the  corresponding  terms  of  the  AWC75  and  IGS75  models. 

Though  a  continuity  in  time  was  achieved  by  the  use  of  secular  change  coeffici¬ 
ents  of  IGRF65  in  the  production  of  IGRF75,  the  new  model  is  not  a  good  represents 
uivw'  of  the  main  field  even  at  the  epoch  1975.0  because  of  the  inherent  errors  of 
the  IGRF65  model.  Accordingly  a  disclaimer  was  adopted  at  the  IAGA  meeting 
stating  that  the  reference  field  is  not  intended  as  a  source  of  compass 
information  for  nautical  and  aeronautical  charts.  Despite  its  poor  represen¬ 
tation  of  the  main  field,  IGRF  models  are  still  very  widely  used  in  the 
reduction  of  geophysical  survey  data. 

AWC70  (Hurwitz  et  al,  1974)  and  AWC75  (Peddle  and  Fabiano,  1976)  models 
were  derived  primarily  for  the  preparation  of  the  American  World  Charts  for 
the  epoch  1970.0  and  1975.0  respectively.  The  main  field  terms  in  both  are 
of  maximum  degree  and  order  12.  The  time  terms  were  basically  derived  from 
observatory  data,  separately  from  the  main  field  analyses,  and  have  maximum 
degree  8. 

The  model  IGS75  was  developed  for  the  preparation  of  British  World  Charts 
(Barraclough  et  al,  1975)  for  epoch  1975.0  and  has  main  field  coefficients 
of  maximun  degree  12.  This  model,  in  addition  to  secular  variation  terms  also 
has  secular  acceleration  terms. 

The  most  recent  model  WC  80  (Barraclough  and  Barker,  1980)  is  the  product 
of  a  joint  effort  by  the  United  States  and  United  Kingdon.  It  is  to  be  used 
to  construct  the  1980  World  Magnetic  Delcination  Charts  published  by  the  Defense 
Mapping  Agency  Hydrographic/Topographic  Center  in  the  U.S.  and  the  British 
Hydrograp.iic  Department  in  the  U.K.  The  main  field  coefficients  are  of  degree 
and  order  12,  and  the  secular  field  coefficients  are  of  degree  and  order  8. 


Models  P0GO2/72  and  GSFC8/73  have  been  derived  primarily  from  the  POGO 
series  satellite  data  and  have  considerably  higher  main  field  and  secular  vari¬ 
ation  coefficient  terms.  The  time  terms  in  these  models  were  determined  simul¬ 
taneously  with  the  main  field  terms  in  the  same  least  squares  analysis. 

Model  MGST3/8Q  (Langel  et  al. ,  1980)  ie  a  preliminary  model  derived  from 
just  two  days  (November  5  and  6,  1979)  of  data  from  recently  launched  satellite, 
MAGSAT.  The  model  is  of  particular  interest,  as  it  is  the  first  time  satelliter 
measured  scalar  and  vector  component  data  of  the  main  field  have  been  included 
in  a  model  construction.  A  better  field  model  from  the  MAGSAT  satellite  data, 
however,  is  expected  to  be  derived  shortly,  as  more  data  are  processed  and  made 
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COMPARISON  OF  FIELD  MODELS 


The  selection  and  use  of  any  particular  field  model  depends  directly  on 
the  requirement  of  a  user,  for  example,  a  user  may  select  a  model  derived  for 
the  epoch  1965.0  if  he  has  to  reduce  survey  data  collected  during  that  time 
period,  whereas  in  another  case,  a  user  may  require  a  model  derived  for  the 
latest  epoch  if  he  has  to  use  it  for  the  purpose  of  navigation.  Difficulties 
arise,  however,  as  there  may  be  several  models  available  for  the  same  epoch 
and  thus  it  becomes  a  tough  choice  to  select  a  proper  model. 

Models  are  primarily  Judged  on  their  prediction  capabilities,  both 
spatially  and  temporally.  Many  papers  have  appeared  in  the  recent  literature 
where  the  performance  of  a  particular  model  and  its  comparison  with  other  models 
have  been  evaluated  in  general  terms.  Such  comparisons  sometimes  may  be  mean¬ 
ingless  or  misleading  as  the  measured  data  against  which  the  model  are  compared 
represent  only  a  limited  sample  and  may  sometimes  even  be  biased  because  of 
local  long  wavelength  anomalies  not  represented  in  the  field  model.  Moreover, 
a  model  judged  superior  to  represent  the  field  in  a  particular  area  may  not  be 
better  or  mav  even  be  worse  in  another  region  of  the  globe.  Nonetheless  results 
of  such  inter-comparisons  provide  valuable  information  about  the  general  utility 
of  a  particular  field  model.  Another  criteria  for  comparison  of  spherical 
harmonic  field  models  is  their  power  spectra.  By  plotting  the  average  degree 
variance  of  the  models,  defined  as 
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one  can  infer  the  depth  of  the  source  mechanism  for  the  various  degree  harmonics. 
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If  these  source  depths  do  not  change  with  time,  then  evaluation  of  the  power  of 
models  at  their  various  epochs  should  yield  the  same  values.  A  model  which  doas 
not  agree  with  the  concensus  in  source  depth  may  be  regarded  with  suspicion. 
Below,  a  brief  summary  of  the  results  of  intercomparison  of  some  of  the  recent 
models  and  their  power  spectra  are  presented. 

The  deficiency  and  limitations  of  the  International  Geomagnetic  Reference 
Field  (IGRF),  a  most  widely  used  model,  are  well  reported  in  the  literature 
(e.g.,  Regan  and  Cain,  1975).  This  model,  as  has  been  discussed  earlier,  is 
a  composite  of  various  models  and  was  chosen  initially  for  epoch  1965.0. 
Barraclough  et  al  (1978)  have  recently  published  and  presented  another  model 
for  epoch  1965.0  which  has  been  derived  from  a  data  set  spanning  the  period  of 
1955  through  1975.  As  the  mean  of  the  data  set  correspond  to  the  epoch  1965.0 
of  the  new  model,  and  more  data  could  not  be  added  to  further  refine  it,  they 
have  referred  to  it  as  a  definitive  model  for  1965.0.  The  authors  have  compared 
their  definitive  model  against  the  IGRF65  model  end  the  other  candidate  models 
from  which  IGRF65  was  selected  and  averaged  and  have  shown  it  to  be  considerably 
superior  in  representing  the  main  field  and  its  secular  variations  around  the 
epoch  of  1965.0.  They  noted  that  though  main  field  coefficients  of  IGRF65  model 
represented  the  main  field  reasonably  well  at  epoch  of  1965.0,  the  poor  secular 
term  coefficients  of  the  model  severely  limited  its  usefulness  beyond  1965.0 
They  further  showed  that  although  IGRF65  model  was  considerably  superior  to 
other  individual  candidate  models,  its  usefulness  would  have  been  considerably 
improved  by  a  more  judicious  selection  end  proper  weighting  of  different  model 
coefficients  than  was  actually  done. 

Mead  (1979)  has  evaluated  performance  of  four  recent  models  namely 
AWC75  (Peddle  and  Fabiano,  1976),  IGS75  (Barraclough  et  al. ,  1975),  1GRF75 
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(IAGA,  1976)  and  POGO  (8/71)  (Langel,  1974),  by  comparing  their  predicted 
values  against  a  set  of  observed  annual  mean  data  obtained  from  obaervatoriua 
around  the  globe.  He  calculated  the  residuals  between  the  observatory  annual 
means  and  the  model  predicted  values  for  the  three  elements  of  the  magnetic 
field  namely  B  (total  field),  I  (inclination)  and  D  (declination)  for  the  years 
1973  through  1976.  His  results  for  the  residuals  AB,  AI  and  AD  cos  I  in  a 
histogram  form  are  reproduced  in  Figures  7  through  9  respectively.  As  has  long 
been  discussed,  the  poorer  quality  of  the  IGRF  model  is  clearly  revealed  here. 
The  other  three  models  can  be  judged  as  equally  well  in  representing  the  main 
field  with  estimated  accuracy  on  the  order  of  150-200  gammas  and  .2°  in  the 
amplitude  and  direction  of  the  field  respectively. 

Coles  (1969)  has  also  made  a  similar  comparative  study  of  the  prediction 
capabilities  of  several  different  recent  models  over  regions  of  Canada.  His 
study  included  seven  different  models  (IGRF65,  TGRF75,  AWC70,  AWC75,  P0GO2/72, 
GSFCS/73  and  1GS75)  and  he  considered  data  from  observatories  as  well  as  from 
aeromagnetic  surveys.  His  results  comparing  model  predicted  data  with  the  ob¬ 
servatory  data  indicate  that  in  general  such  predictions  are  in  better  agreement 
with  the  portion  of  observatory  data  corresponding  to  the  period  of  main  data 
set  used  in  the  production  of  the  field  model.  However,  problems  develop  out¬ 
side  this  time  range  mainly  because  of  the  non-linearity  of  the  rate  of  secular 
change.  In  applying  field  models  to  the  reduction  of  airborne  survey  data,  his 
study  clearly  pointed  out  the  advantage  of  models  having  higher  order  harmonic 
terms  as  they  afford  a  much  better  spatial  resolution  needed  in  such  regional- 
residual  separation. 

Langel  et  al.  (I960)  has  recently  proposed  a  new  model  (MGST  3/80)  based 
on  2-daya  of  MAGSAT  satellite  data.  They  have  also  used  this  set  of  data  to 
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compare  the  prediction  capabilities  of  other  recent  models  at  the  MAGSAT  altitude 
and  epoch.  Table  3  summarily  presents  the  results  of  their  comparative  study. 

The  results  clearly  point  out  that  POGQ2/72  model  which  was  based  only  on  scalar 
measurements  made  by  earlier  POGO  satellites,  is  the  best  model  in  correctly 
predicting  the  scalar  field  measured  by  MAGSAT.  However,  AWC75  and  IGS75 
appear  to  better  describe  the  current  vector  fields. 

Comparison  of  the  power  spectra  of  seven  spherical  harmonic  modelB  was 
carried  out  following  the  procedure  of  Cain  (1975).  The  power  of  a  model  was 
evaluated  using  equation  (5)  and  the  epoch  value  of  the  model  coefficients. 

The  power  spectra  values  for  the  P0CO6/74,  P0G02/72,  AWC75,  IGRF75,  GSFC12/66, 
IGRF10/68,  and  the  WC80  models  are  displayed  in  Figure  10,  along  with  an 
approximate  power  law  curve  fit  adapted  from  Cain  (1975).  The  different  slopes 
of  the  straight-line  segments  of  this  curve  are  interpreted  by  Cain  as  represen¬ 
tative  of  the  core  (n<8) ,  crust  (n>13)  and  mantle  transition  (&<n<13)  source 
region. 

All  of  the  models  yield  virtually  the  same  powers  for  n<8,  and  even  for 
n-10,  no  significant  disagreements  exist.  However,  for  n*ll  and  12,  two  models 
WC80  and  FOG02/72  yield  significantly  higher  power  than  P0G06/74  and  AWC75. 

This  may  be  due  to  the  effects  of  aliasing,  or  truncation  error,  since  the 
higher  degree  P0G06/74  model  has  much  lower  power  at  n*12.  Tlie  sudden  increase 
in  power  at  n*13  for  P0G02/72  is  especially  suspicious.  Since  the  WC80  and 
POC06/71  models  both  have  epochs  which  are  many  years  away  from  those  of  the 
P0G06/74  or  AWC75  models,  it  is  possible  that  the  observed  differences  in 
power  are  somehow  due  to  secular  changes  in  the  dapth  of  the  source  mechanisms. 
This  question  may  be  resol’d,  by  study  of  additional  models  having  common  epochs. 
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normalized  power  in  spherical  harmonics  of  internal  field  models 
evaluated  at  the  appropriate  epoch.  The  power  law  curve  fit  is 
taken  from  Cain  (1975). 


SUMMARY  AND  CONCLUSIONS 


The  Earth's  magnetic  field  is  complex  and  dynamically  varying  vith 
approximately  99X  of  the  field  at  the  Earth's  surface  originating  from  internal 
sources  and  the  remainder  from  sources  outBide  the  Earth,  mainly  in  the  ionosphere 
and  magnetosphere.  That  part  originating  from  internal  sources  is  called  the 
internal  field  or  the  main  field  and  that  originating  from  external  sources 
referred  to  as  the  external  field.  The  main  field  varies  relatively  smoothly 
over  the  Earth's  surface  and  is  approximately  dipolar  in  nature.  The  main 
field,  however,  is  not  permanent  in  time  but  rather  varies  relatively  slowly 
over  period  of  years  and  such  variation  in  the  main  field  is  termed  secular 
variation. 

From  measurements  of  the  field  over  the  Earth's  surface,  mathematical 
models  of  the  main  field  and  its  secular  variation  can  be  derived  which  describe 
the  measured  field.  Such  a  mathematical  representation  of  the  global  main  field 
is  usually  done  in  terms  of  a  truncated  spherical  harmonic  series  and  is  called 
a  geomagnetic  field  model  or  simply  field  model.  Such  models  are  essentially 
data  extrapolation  functions,  i.e.,  from  a  set  of  measured  global  data  they  are 
used  in  to  extrapolate  field  values  both  spatially  and  temporally.  This  analyt¬ 
ical  representation  of  the  global  main  field  makes  field  models  ideally  suited 
for  data  processing  and  are  used  in  a  variety  of  applications  including  preparation 
of  world  magnetic  charts  for  navigational  purposes,  removal  of  regional  trends 
from  magnetic  surveys  to  isolate  local  crustal  anomalies,  and  in  numerous  space 
science  studies. 

Comparison  and  evaluation  of  recent  field  models  indicate  that  the  best 
field  models  could  describe  the  field  vith  overall  accuracies  batter  than  approx¬ 
imately  150-200  gammas  in  magnitude  and  .2*  in  direction  of  the  field.  Effects 


of  crustal  sources  as  well  as  the  external  sources  are  not  Included  in  the 
field  models  due  to  their  limited  spatial  and  temporal  frequencies.  This  may 
be  a  limitation  of  current  field  models  and  continued  improvement  in  accurate 
representation  of  the  Earth's  field  would  require  inclusion  of  these  two 
sources,  namely,  the  crustal  sources  and  the  external  sources  in  the  field 
models.  Inclusion  of  crustal  anomalies  in  the  model,  however,  la  neither  war¬ 
ranted  nor  desired  for  geophysical  applications  but  may  be  highly  desirable  for 
other  applications  such  as  navigation  and  space  science  studies  Inclusion  of 

such  a  signal  would  require  carrying  out  the  spherical  harmonic  series  analysis 
to  a  much  higher  degree  and  order  than  done  at  present.  This  is  a  formidable 
task,  not  only  because  of  the  immensely  Increased  computer  time  and  storage 
requirement  but  also  because  the  necessary  detailed  global  spatial  data  are 
not  available.  A  better  solution  might  be  to  add  local  function  representa¬ 
tions  of  the  crustal  anomalies  of  a  particular  region  to  the  global  field 
model,  which  then  would  allow  a  much  better  detailed  definition  of  the  field 
over  that  particular  region. 

Another  problem  in  accurate  representation  of  the  geomagnetic  field  by 
field  models  comes  from  the  fact  that  the  Earth's  field  is  in  a  continual  state 
of  change.  Such  secular  change  must  be  monitored  and  retained  in  the  model,  if 
model  prediction  Is  to  be  valid  over  future  years.  Lack  of  enough  data  has  gen¬ 
erally  prevented  accurate  global  modeling  of  secular  changes  and  only  first 
order  linear  time  variation  terms  over  the  period  of  the  data  base  are  calculated 
and  retained  in  current  models.  Few  models  containing  second  order  time  terms 
(secular  acceleration)  have  been  calculated.  Such  second  degree  time  terms  may 
in  certain  instances  improve  the  description,  but  can  also  worsen  the  predictive 
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frequent  updating  of  field  models  is  necessary  to  account  for  the  secular 
changes  in  the  Earth's  main  field. 

Mapping  of  Earth’s  magnetic  field  by  satellites  in  the  past  decade  have 
considerably  aided  in  better  representation  of  the  field  as  they  provided  es¬ 
sentially  global  coverage  over  a  very  short  period  of  time.  In  the  past,  such 
data  were  limited  to  the  total  field  measurements  only  and  as  Stern  and  Bredekaap 
(1974)  have  noted  the  models  based  on  these  measurements  may  be  very  pre¬ 
cise  in  representing  total  field  values  at  satellite  altitudes,  but  they  may  be 
deficient  in  determining  accurate  vector  values.  Such  deficiencies,  however, 
are  expected  to  be  removed  from  the  analysis  of  MAGSAT  satellite  data,  which 
is  to  provide  vector  field  measurements  as  well  as  the  total  field  values. 

Despite  such  satellite  data,  the  role  of  airborne  survey  data  like  Project 
MAGNET  data  in  defining  the  Earth's  field  cannot  be  denied  as  the  two  types  of 
data  are  essentially  complimentary  to  each  other  and  analysis  of  combined  data 
sets  result  in  a  much  better  model  of  the  Earth's  field. 
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